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- - Abstract 

, We consider J /ip meson photoproduction on protons at high energies at the leading 

I order in using conventional parton model and /cy-factorization approach of QCD. It 

pH ■ is shown that in the both cases the colour singlet mechanism gives correct description 

a , for experimental data from HERA for the total cross section and for the J/ip meson 

z-spectrum at realistic values of a c-quark mass and meson wave function at the origin 
^'(0). At the same time our predictions for pf-spectrum of J/ip meson and for py 
dependence of the spin parameter a obtained in A:T-factorization approach are very 
\ different from the results obtained in conventional parton model. Such a way the 

experimental study of a polarized J/'i/' meson production at the large pj should be a 
CN ■ direct test of BFKL gluons. 

o 



X 



^ : 1 Introduction 
o 

^ ' R is well known that in the processes of J/ip meson photoproduction on protons at high 

energies the photon-gluon fusion partonic subprocess dominates |jl|. In the framework of 
general factorization approach of QCD the J/ip photoproduction cross section depends on 
r-| ! the gluon distribution function in a proton, the hard amplitude of cc-pair production as 

well as the mechanism of a creation colorless final state with quantum numbers of the J/ip 
meson. Such a way, we suppose that the soft interactions in the initial state are described by 
introducing a gluon distribution function, the hard partonic amplitude is calculated using 
perturbative theory of QCD at order in ^^(/i^), where /x ~ mc, and the soft process of the 
cc-pair transition into J/ip meson is described in nonrelativistic approximation using series 
in the small parameters as and v (relative velocity of the quarks in J/ip meson). As it is 
statement in nonrelativistic QCD (NRQCD) 0, there are colour singlet mechanism, where 
cc-pair is hardly produced in the colour singlet state, and colour octet mechanism, where 
cc-pair is produced in the color octet state and at long distance it transforms into final colour 
singlet state in soft process. However, as it was shown in paper P], the data from HERA 
collider in the wide region of pt and z may be described well in the framework of the colour 
singlet model and not colour octet contribution is needed. Based on the above mentioned 
result we will take into account in our analysis only the colour singlet model contribution 
in the J/ip meson photoproduction [jl|. We consider a role of a proton gluon distribution 
function in the J/ ip photoproduction in the framework of the conventional parton model as 
well as in the framework of the fcy-factorization approach ||^. Last one is based on BFKL 
evolution equation |]^ , which takes into account large terms proportional to log 1/x and 
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log/i^/Ag^^ opposite to DGLAP evolution equation 0, where only large logarithmic terms 
~ log/x^/Ag^^ are taken into account. In the process of the J/ip meson photoproduction 
one has x ~ /i^/s and //^ ~ M^, where M is the J/ip meson mass, s is the square of a total 
energy of colliding particles in center of mass frame. 

2 The cross section for J/ipX in /cy-factorization 

approach 

Nowadays, there are two approaches in calculating of J / ip meson or heavy quark production 
cross sections at high energies. In the conventional parton model P it is suggested that 
hadronic cross section o"(7p — J /ipX^ s) and relevant partonic cross section (y{'-)g — > J /ipg^ s) 
are connected as follows: 

aP'\^p^ J/,IjX,s)= ldxa{jg^J/ijg,s)G{x,fi^), (1) 

where s = xs, G{x,ii^) is the collinear gluon distribution function in a proton, x is the 
fraction of a proton momentum, /x^ is the typical scale of a hard process. The fi^ dependence 
of the gluon distribution is described by DGLAP evolution equation 0. In the 

region of very high s one has a; << 1. This fact leads to BFKL evolution equation for the 
unintegrated gluon distribution function $(x,qj,yU^), where qy is the gluon virtuality. The 
unintegrated gluon distribution function can be related to the conventional gluon distribution 
by 

a;G(x,/i^)= / qy,;U^)(iq?-. (2) 

Jo 

The gluon 4-momentum is presented as follows 

q = xpN + qj, 

where = (0,qT,0), pn = {En,0,0, \pn\) and q^ = q\ = — qr^- So called BFKL gluon 
is off mass-shell and it has polarization vector along its transverse momentum such as e'^ = 
Qt/\^t\- In fcy-factorization approach hadronic and partonic cross sections are related by 
the following condition: 

a{^p^J/^X) = [ — r [ ^^x,q\,fi')a{^g* ^ J/^g,s,q\), (3) 
J X Jo J ZTT 

where d'{'yg* J/'ipo) is J/'ip meson photoproduction on BFKL gluon, is azimuthal angle 
in the transverse XY plane between vector qy and fixed OX axis. 



3 Unintegrated gluon distribution function 

At the present time an exact form of the unintegrated gluon distribution $(x,qy,;U^) is 
unknown yet because the relevant experimental analysis of the experimental data has never 
been carried out. There are several theoretical approximations for qy,/i^), which are 
based on solving BFKL evolution equation [1^. In the region of very small x < 0.01 

and moderate q^ (~ 10 GeV^), which is relevant for J/^f photoproduction at HERA, all 
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parameterizations |T0|, |T^ are approximately coincide (see, for example, discussions in 



H\ and |T^). For our purposes, we will use the unintegrated gluon distribution, which was 
obtained in the paper [|I4|. The proposed method lies upon a straightforward perturbative 
solution of the BFKL equation where the coUinear gluon density G{x,iJ?) is used as the 
boundary condition in the integral form (2). Technically, the unintegrated gluon distribution 
is calculated as a convolution of coUinear gluon distribution G'(x,yU^) with universal weight 
factors: 

$(x,q^,/i2) = r6;(r^,q^,/x2)-G(-,/i2)rfr/, (4) 
Jx n n 



where 




a 



r/q 



as In(-) ln(^; 
r] jj,'^ 



, at qj < /i^ 



at qy > /i^. 



(5) 
(6) 



Jo and Jo stand for Bessel functions (of real and imaginary arguments, respectively), and Us = 
— . As a input function G(x, jj?) in (5,6) we use the standard GRV parameterization [|15|. To 

test the method of calculation for $(x, q^, /i^) we compare input coUinear gluon distribution 

(^(x, /i^) and coUinear gluon distribution ^(x, /x^), which is obtained using formula (2) from 

the unintegrated distribution $(x, qy,/i^) after integration over q^. In the Figure 1 the 

result of our calculation for the ratio i?i(x, /i^) = G{x, fi'^)/G{x, fi"^) is shown. It is visible, 

that at X < 0.1 the ratio Ri does not differ from 1 more than 1-2 %. Note, that in a process 

M 

of the J/ip photoproduction on protons x ~ —= (M is the mass of J/ip meson), and one has 
X ~ 0.03 at = 100 GeV. 



4 Amplitude for J/i/^g process 

There are six Feynman diagrams (Figure 2) which describere partonic process 'jg J/'4'9 

at the leading order in and a. In the framework of the colour singlet model and nonrel- 

ativistic approximation the production of J/ip meson is considered as the production of a 

quark- ant iquark system in the colour singlet state with orbital momentum L = and spin 

momentum S = 1. The binding energy and relative momentum of quarks in the J/ip are 

P 

neglected. Such a way M = 2mc and Pc = Pc = -, where p is 4-momentum of the J/V", Pc 
and Pc are 4-momenta of quark and antiquark. Taking into account the formalism of the 



projection operator |jT6| the amplitude of the process 7(7 J/4'9 be obtained from the 



amplitude of the process 'jg ccg after replacement: 

V\p,)W{pc) - ^^iP)iP + M)^, (7) 

where i{p) = s^{p)'~^^, ^^l{p) is 4-vector of the J/il) polarization, —= is the color factor, \E'(0) 

V3 

is the nonrelativistic meson wave function at the origin. The matrix elements of the process 
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^g* — > J/i^g may be presented as follows: 

M, = KC^'e^{k,)el{q)e'p{k2)eMMf'\ 



(8) 



m: 



Tr 
Tr 
Tr 
Tr 
Tr 
Tr 



- , nj\ Pc-ki + rric -pc - h + nic a 

7 p + M 7 — _^v.\2 2^ 

[Pc — kiY — [pc + ^2) — fn^ 

7^(p + M)7'^-^'= + ^^ + '"'= g-Pc- + m, 



r7 



{pc + — ml {q — pc)^ — ml 



r 



i^f- , a Pc-h + mc 8 q-Pc + m^ 



(Pc — — ml {q — PcY ~ f^t 



Y{p + M)Y 



^ Pc-q + mc -pc -k2 + mc 
7"(p + M)7^ ' ' -^'^ "'^ ' 



7'^(p + M)7^ 



(Pc - g)^ - T^ic ' (^"2 + Pc)^ - m\ 

Pc + h + mc „ ki-Pc + mc 
7 7 

(Pc + ^2)^ — ""^c (^1 ~ Pc)^ — ''"c 

Pc-q + mc ki-pc + mc 



{pc-qY-ml {ki-pcY-m: 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 



where ki is the 4-momentum of the photon, q is 4-momentum of the initial gluon, k2 is the 
4-momentum of the final gluon, 



K = Cccgl 
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C"^^ = -^Tr [T'^T^], ee=| e = x/i^, = x/4^. 



The summation on photon, J/ip meson and final gluon polarizations is carried out by co- 
variant formulae: 



Y2sa{ki)e/3{ki) = -ga(3, 

spin 

spin 



spin 



M2 ■ 



(15) 
(16) 
(17) 



In case of the initial BFKL gluon we use the following prescription 



Y,Si,{q)e^{q) 



spin 



(18) 



For studing J/ip polarized photoproduction we introduce the 4- vector of the longitudinal 
polarization as follows 



pi" Mp^j^ 



(19) 



M (ppn) ■ 

In the limit of s >> the polarization 4- vector satisfies usual conditions (elSl) = — 1, 
(slp) = 0. 
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Traditionally for a description of quarkonium photoproduction processes the invariant 
variable z = (pPn)/ (kipiy) is used. In the rest frame of the proton one has z = E^/E^. In 
the fcy-factorization approach the differential on pj- and z cross section of the J/ip photo- 
production may be written as follows 



1 — v/ ^ t^qT'^'(a;,qT,/^ )T7ry^- (20) 

1 — z) Jo In Jo loTT xs r 



dp\dz z{) 

The numerical calculation is performed in the photon and proton center of mass frame where 

PN = ^(1,0,0,1), h = ^(1,0,0,-1), 

p = (^,Pt,P||), g = (— x,qT,— x). 

Here we take into account that the J/ip momentum p lies in XZ plane and (qtPt) = 
IptHqtI cos(0). The analytical calculation of the |Mp is performed with help of REDUCE 
package and results are saved in the FORTRAN codes as a function of s = {ki + g)^, 
t = {p — kiY, u = {p — q)^, py, Qy and cos(0). We directly have tested that 

hm r^W? = W^, (21) 

iu 



where py = — in the |Mp and |MpA/P is the square of the amplitude in the conventional 



s 



parton model 0]. In the limit of qj = from formula (20) it is easy to find the differential 
cross section in the parton model too: 



dp\dz \Qt[{xsYz{\ — z) 

However, making calculations in the parton model we use formula (20), where integration 
over qy and is performed numerically, instead of (22). This method fixes common normal- 
ization factor for both approaches and gives direct opportunity to study effects connected 
with virtuality of the initial BFKL gluon in the partonic amplitude. 



5 Results and discussion 

After we fixed selection of the gluon distribution function G{x, /i^) there are two parameters 
only, which values determine the common normalization factor of the cross section under 
consideration: \l/(0) and rric. The value of the J/ip meson wave function at the origin may 
be calculated in a potential model or obtained from experimental well known decay width 
V{J/ip — > fj,~^fj,~). In our calculation we used the following choice |\l/(0)p = 0.0876 GeV^ 

as the same as in Ref. 0. Concerning a charmed quark mass, the situation is not clear 

M 

up to the end. From one hand, in the nonrelativistic approximation one has rric = — , but 

there are many evidences to take smaller value of a c-quark mass in the amplitude of a hard 
process, for example rric =1.4 GeV. Taking into consideration above mentioned we perform 
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calculations as at rric = 1.4 GeV as well as at = 1.55 GeV. The cinematic region under 
consideration is determined by the following conditions: 0.4 < z < 0.9 and p\ > 1 GeV^, 
which correspond HI and ZEUS Collaborations data We assume that the contribution of 
the colour octet mechanism is large at the z > 0.9 only. In the region of the small values of 
z < 0.2 the contribution of the resolved photon processes [jl^ as well as the charm excitation 



processes |18] may be large too. All of these contributions are not in our consideration. 



i?2(S,q^T) = ^^, (23) 



Figure 3 shows the ratio 

It 

a(s,0) 

as a function of at the different s. As well as it was necessary to expect ratio (23) decrease 
with growth of and the faster than less the value of s. 

Figures 4-6 show our results which were obtained as in the conventional parton model 
as well as in the /cy-factorization approach at two values of a charmed quark mass. The 
dependence of the results on selection of a hard scale parameter /i is much less than the 
dependence on selection of a c-quark mass. We put /i^ = + py in a gluon distribution 
function and in a running constant ^^(/i^). 

Figure 4 shows the dependence of the total J/ip photoproduction cross section on y/s. 
It is visible that the difference between the parton model prediction and the result of kj- 
factorization approach is much less than between the results obtained at different values of 
c-quark mass in the both models. At the rric = 1.55 GeV the obtained cross sections in 1.5 
-2 times are less than experimental data [Q, but at the rric = 1.4 GeV the theoretical curves 
lie even a shade higher of the experimental points. 

There are not contradictions between the theoretical predictions and data for the z- 
spectrum of the J/'0 mesons. Figure 5 shows that the experimental points lie inside the 
theoretical corridor as in the parton model as in the fcy-factorization approach. 

The count of a transverse momentum of the BFKL gluons in the /cj-factorization ap- 
proach results in a flattening of the py-spectrum of the J/ ip as contrasted to by predictions 
of parton model. For the first time on this effect was indicated in the article |]T9|, and later 



in the articles [T^, |2^. Figure 6 shows the result of our calculation for the py-spectrum 
of the J/ip mesons. Using the /cy-factorization approach we have obtained the more hard 
PT-spectrum of the J/tp than it have been predicted in the LO parton model. It is visible 
that at large values of pj only the /cy-factorization approach gives correct description of the 
data . However, it is impossible to consider this visible effect as the direct indication on a 
nontrivial developments of the small-x physics. In the article |^ was shown that calculation 
in the NLO approximation gives more hard py-spectrum of the J/ip meson too, which one 
will be agreed the data at the large pj- 

As it was mentioned above, the main difference between the fcy-factorization approach 
and the conventional parton model is nontrivial polarization of BFKL gluon. It is obviously, 
that such spin condition of the initial gluon should result in observed spin effects during 
birth of the polarized J/tp meson. We have performed calculations for the spin parameter a, 
as a function z or p^ in the conventional parton model and in the /cy-factorization approach 

datot g dax^ 

dz dz 
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dcrtot g c^cl 

"(Pt) = tL , ^ (25) 

Here atot = ctl + ax is the total J/ip production cross section, cr^ is the production cross 
section for the longitudinal polarized J/ip mesons, ax is the production cross section for 
the transverse polarized J/ip mesons. The parameter a controls the angle distribution for 
leptons in the decay J/tp —>■ in the J/ip meson rest frame: 

l + acos(0). (26) 



dcos{6) 



Figure 7 shows the parameter a{z), which is calculated in the parton model (curve 2) and 
in the /cy-factorization approach (curve 1). We see that both curves lie near zero aX z < 0.8 
and increase sX z > 0.8. The large difference between predictions is visible only aX z > 0.9 
where our consideration is not adequate. Let's remark, that parameter is gentle depends on 
mass of a charmed quark and we demonstrate here only outcomes obtained at the tjIc = 1.55. 

For the parameter a{pT) we have found strongly opposite predictions in the parton model 
and in the fcy-factorization approach, as it is visible in Figure 8. The parton model predicts 
that J/ if) mesons should have transverse polarizations at the large pt (^(pt) = 0.6 at p\ = 20 
GeV^), but fcy-factorization approach predicts that J/ip mesons should be longitudinally 
polarized ( ^(pt) = —0.4 at pj = 20 GeV^). Nowadays, a result of NLO parton model 
calculation in the case of the polarized J/ip meson photoproduction is unknown. It should 
be an interesting subject of future investigations. If the count of the NLO corrections will 
not change predictions of the LO parton model for a{pj), the experimental measurement of 
this spin effect will be a direct signal about BFKL gluon dynamics. 

Nowadays, the experimental data on J/\E' polarization in photoproduction at large pj 
are absent. However there are similar data from CDF Collaboration I^T], where J/ip and 
ip' pT-spectra and polarizations have been measured. Opposite the case of J/ip photopro- 
duction, the hadroproduction dada needs to take into account the large colour-octet contri- 
bution in order to explain J/ip and ip' production at Tevatron in the conventional coUinear 
parton model. The relative weight of colour-octet contribution may be smaller if we use 



fcy-factorization approach, as it was shown recently in ||22| , |23| , p^ , |25| . The predicted using 
collinear parton model transverse polarization of J/ip at large pj is not supported by the 
CDF data, which can be roughly explain in /cy-factorization approach [^. In conclusion, 
the number of theoretical uncertainties in the case of J/ip meson hadroproduction is much 
more than in the case of photoproduction and they need more complicated investigation, 
that is why the future experimental analysis of J/ip photoproduction at THERA will be 
clean check of the collinear parton model and fcy-factorization approach. 

This work has been supported in part by the programme "Universities of Russia - Basic 
Researches" under Grant 02.01.003. 
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Figure 1: Diagrams used for description partonic process 7 + 5'— > J/ip + g. 
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Figure 2: Ratio -Ri(a;,/i2) at /x^ = 10 GeV^ (curve 1) and /x^ = 100 GeV^ (curve 2). 
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a , nb 




Figure 4: The J/ip photoproduction cross section as a function of ^/s at 0.4 < z < 0.9 and 
p\ > 1 GeV^. The data from [4], the curves 1 and 2 are obtained in the /cy-factorization 
approach at rric = 1.4 and rric = 1.55 GeV, the curves 3 and 4 are obtained in the parton 
model at rric— 1.4 GeV and rric — 1.55 GeV. 



11 



d a/dz, nb 




.0 0.2 0.4 0.6 0.8 1.0 



z 

Figure 5: The J /ijj spectrum on z at ^/s = 100 GeV and p\ > 1 GeV^. The data from [4], 
the curves are the same as in Figure 4. 
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da/dp2 ,nb / GeV 
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Figure 6: The J/ip spectrum on p\ at \/s = 100 GeV and 0.4 < z < 0.9. The data from [4], 
the curves are the same as in Figure 4. 
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Figure 7: The parameter a as a function of z at \fs = 100 GeV, p\ > 1 GeV^ and rric = 1.55 
GeV. The curve 1 is the result obtained in the /cT-factorization approach, the curve 2 is the 
result obtained in the parton model. 



14 



a (p2) 




